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ABSTRACT 

Radio polarimetry at decimetre wavelengths is the principal source of information on the Galactic magnetic 
field. The diffuse polarized emission is strongly influenced by Faraday rotation in the magneto-ionic medium 
and rotation measure is the prime quantity of interest, implying that all Stokes parameters must be measured 
over wide frequency bands with many frequency channels. The DRAO 26-m Telescope has been equipped 
with a wideband feed, a polarization transducer to deliver both hands of circular polarization, and a receiver, all 
operating from 1277 to 1762 MHz. Half -power beamwidth is between 40 and 30 arcminutes. A digital FPGA 
spectrometer, based on commercially available components, produces all Stokes parameters in 2048 frequency 
channels over a 485-MHz bandwidth. Signals are digitized to 8 bits and a Fast Fourier Transform is applied to 
each data stream. Stokes parameters are then generated in each frequency channel. This instrument is in use 
at DRAO for a Northern sky polarization survey. Observations consist of scans up and down the Meridian at a 
drive rate of ^0.9° per minute to give complete coverage of the sky between declinations —30° and 90°. This 
paper presents a complete description of the receiver and data acquisition system. Only a small fraction of the 
frequency band of operation is allocated for radio astronomy, and about 20% of the data are lost to interference. 
The first 8% of data from the survey are used for a proof-of-concept study, which has led to the first application 
of Rotation Measure Synthesis to the diffuse Galactic emission obtained with a single-antenna telescope. We 
find rotation measure values for the diffuse emission as high as ~ ±100 rad m~ 2 , much higher than recorded 
in earlier work. 

Subject headings: Instrumentation: polarimeters - methods: data analysis - ISM: magnetic fields - polarization 
- surveys 



1. INTRODUCTION 

All extensive surveys of the Galactic polarized emission 
to date have measured polarization properties of the emis- 
sion at one or a few frequencies over narrow bandwidths. 
Recent examples at decimetre wavel engths are the DRAO 
Low-Resolution Polarization Survey ( Woll eben et al.ll2006l) . 
the Argentinian polarization survey dTestori et al. 2008), and 
the Effelsberg Medium Latitude Survey (iReich et aL 2004). 
However, it is widely accepted that the appearance of the po- 
larized sky at these wavelengths depends less on the signal 
properties at the point of emission and more on Faraday ro- 
tation along the propagation path; that appears to be true for 
the Milky Way emission at all frequencies up to about 5 GHz. 
If polarization surveys do not contain information on rotation 
measure (RM), the quantity that is significant for physical un- 
derstanding of the interstellar medium (ISM), their value is 
limited. 

Recent years have seen the development and application of 
wide-band feeds and receiver systems. Advances in signal- 
processing technology have led to construction of digital full- 
polarization correlators, often based on Field-Programmable 
Gate Arrays (FPGAs) or similar devices. In combination with 
commercial analog-to-digital converters (ADCs) with high 
sampling rates, often exceeding 1 GHz, FPGAs can be used 



to realize digital polarimeters with large input bandwidths and 
high dynamic range. Wideband spectro-polarimetric data can 
now be easily acquired. 

Methods for the an alysis o f such data were articulated over 
four decades ago by iBurnl d 19661) but the technical limita- 
tions of the time made the observations and the ir analysis ex- 
tremely difficult. Brentje ns & de Bruynl (120051) have recently 
developed Rotation Measure Synthesis (RM-Synthesis), and 
it has been shown to be a p owerful tool for the analysis 
of wide-band polarizatio n data dde B ruvn & Brenti ens 1 120051 : 
ISchnitzeler etaT] l2007allR 120091: iBernardi et al] 120091) . The 
application of this technique requires correctly calibrated, 
multi-frequency polarization data and has to date been exclu- 
sively applied to observations obtained with synthesis tele- 
scopes. 

The output from multi-frequency polarization observations 
is a data cube of images at many wavelengths, A. Applica- 
tion of RM-Synthesis transforms this data cube into an RM- 
Synthesis cube. The first and second dimensions of an RM- 
Synthesis cube are the coordinates on the sky, which are R.A. 
and DEC for the data here presented. The third dimension 
is Faraday depth, <f> [rad m~ 2 ], which is often used in RM- 
Synthesis to replace RM: 
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where By [/iG] is the magnetic field component parallel to the 
line-of-sight, n e [cm -3 ] the electron density, and I [pc] the dis- 
tance along the line-of-sight. Choosing a value for cf) the ob- 
served polarization vector in each frequency channel is rotated 
by </>A 2 and the rotated vectors are coherently added to obtain 
an image of polarized intensity at that value of Faraday depth. 
It is important to note that Faraday depth is usually not related 
to physical distance. 
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Significant wideband spectro-polarimetric studies have 
been made with aperture-synthesis telescopes (see references 
above), but absence of information on zero-levels in Stokes U 
and Q causes non-linear effects in polarization angle and po- 
larized intensity, making interpretation always difficult, and 
sometimes impossible. Incorrect zero levels can turn polar- 
ized objects into apparently unpolarized ones, and can ren- 
der RM values for the diffuse Galactic emission meaning- 
less. If the region mapped is of significant extent aperture- 
synthesis observations must be complemented with absolutely 
calibrated single-antenna data to provide correct zero levels 
and information on the largest structures. 

Taking advantage of the technical developments described 
above, a new generation of polarization surveys with single 
antennas is underway, bringing the field into the era of 3-D po- 
larimetry. Multi-frequency polarization surveys with single- 
antenna telescopes include GALFA CTS with the Arecibo 
telescope from 1225 to 1525 MHz dGibsonl 120051) . STAPS 
with the Parkes telescope from 1296 to 1804 MHz (PI M. 
Haverkorn), and S-P ASS also with the P arkes telescope from 
2188 to 2412 MHz dCarretti et al.ll2008l) . 

The work described here is part of a new initiative, the 
Global Magneto-Ionic Medium Survey (GMIMS), a major 
project to survey the entire sky from the northern and south- 
ern hemispheres covering 300 MHz to 1.8 GHz with many 
thousands of frequen cy channels using large single antennas 
dWolleben et alJ l2008h . The focus of this paper is on the 
technical realization of spectro-polarimetry using the 26-m 
Telescope at the Dominion Radio Astrophysical Observatory 
(DRAO). The telescope has been equipped with a wide-band 
receiver and an FPGA-based FFT polarimeter. The data gath- 
ered to date have led to the first application of RM-synthesis to 
single-antenna polarization data, and we present those results 
and give a general interpretation. The data presented form the 
proof-of-concept study for the survey now in progress. 

2. TELESCOPE, FEED AND RECEIVER 
2. 1 . Reflector and feed 

The axi-symmetric reflector, of diameter 25.6-m, is equato- 
rialy mounted. The mesh surface has an accuracy of ~0.5 cm 
rms. The receiver is mounted at the prime focus on a fibre- 
glass tripod. The telescope can observe the sky between —30° 
and +90° in an almost fully automated fashion. 

The feed is scaled from the design of Wohllebe n et al.1 
(1972), which was developed for the Effelsberg 100-m Tele- 
scope. The DRAO 26-m Telescope has the same ratio of focal 
length to diameter, f/d — 0.3, and the feed provides excellent 
performance at its design frequency, 1420 MHz, where the 
aperture efficiency achieved is ~ 0.53 and the beam is circu- 
lar with half -power width of 36'. 

Fee d performanc e has been computed across the operating 
band (Smegal 2005). Feed radiation patterns are highly circu- 
lar (with equal half-power widths in the E- and //-planes) at 
1400 MHz. At lower frequencies the //-plane pattern is wider 
than the E, and conversely for frequencies above 1400 MHz. 
This leads to increased instrumental polarization at frequen- 
cies away from 1400 MHz. Figure 1 shows the beamwidth 
of the telescope across the entire frequency band. The small 
variations in beamwidth are real (they are highly repeatable) 
but their cause is unknown. 

2.2. Polarizer 

The receiving system is designed to accept all the power 
from an incoming signal and to produce outputs that describe 
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FIG. 1 .— The half-power beamwidth (HPBW) of the DRAO 26-m tele- 
scope for RR (continuous) and LL (dashed) correlation products determined 
by fitting Gaussian brightness distributions to a few observations of our pri- 
mary and secondary calibrators. 

its polarization state. Two antennas of orthogonal polariza- 
tion, whether linearly or circularly polarized, will together 
collect all the incoming power. Given the a priori knowledge 
that the linearly polarized component of sy nchrotron emission 
dominates any circularly polarized content ( Pach olczvkl 1 9701) 
we designed the receiving system to measure small linearly 
polarized components in largely unpolarized signals. The two 
channels of the receiver accept right-hand circular polariza- 
tion (R) and left-hand circular polarization (L). Stokes pa- 
rameters Q and U that describe linear polarization can then 
be derived by correlation techniques. 

Consider the incoming signals as the sum of two circularly 
polarized field components 



E r = E R e' 



and 



E, = E L e iwt+S . 



(2) 



(3) 



In terms of these components the desired outputs from the 
polarimeter are 



/=i(£ L 2 +£ x 2 



Q = ElErCO&5, and 



U = ElErSiiiS. 



(4) 
(5) 
(6) 



V = j(/?l 2 —Er 2 ) is small and we will ignore it. 

The alternative scheme, receiving two linear polarizations, 
leads to the derivation of Q as the (small) difference between 
two large quantities, and is therefore quite susceptible to in- 
strumental effects such as gain variations between receiver 
channels. 

The operation of the feed as a collector of circularly po- 
larized signals can be more easily understood by consider- 
ing it as a transmitter and showing that it radiates circularly 
polarized signals. An orthomode transducer combines two 
monochromatic signals A and B of equal amplitude from 
coaxial inputs, launches them as separate signals in rectan- 
gular waveguide and then combines them in a single square 
waveguide where they travel in orthogonal (non-interacting) 
TEio and TEqi modes. A power splitter in square waveguide 
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converts A and B into \{A+B) and \{A — B). A differen- 
tial phase shifter, also built in square waveguide, provides a 
90° relative phase shift between these signals and so generates 
outputs j(A + iB) and \{B — iA), which are circularly polar- 
ized components. A smooth transition from square to circular 
waveguide transfers these two outputs to the feed, from which 
they will emerge as circularly polarized signals in free space. 

By reciprocity this device will collect the two circularly po- 
larized components E/ and E r described in equations 2 and 3 
from an incoming wave and will deliver them to the coaxial 
inputs of the two receiver channels. In the following we de- 
note E\ simply as L and E r as R, and recognize that they are 
broadband signals. Powers will be denoted as: RR — E r E* , 
LL = EiEf, RL = Re{E r Ef), and LR = Im(E r Ef). 

The return loss of the polarizer at the receiver inputs is 20 
dB: only 1 % of the power in a signal injected into the polarizer 
at either port returns to that port. The isolation between the 
two ports is better than 30 dB (0.1% of signal power). The 
dissipative loss is less than 2% (0.1 dB) and the axial ratio 
of the device (measured by transmitting from it and analysing 
the radiated signal) is of the order of ±10%. 

2.2.1. Phase Shifter, Power Splitter, and Orthomode 
Transducer 

The phase shifter is built in square waveguide (12.7 x 
12.7 cm); its overall length is ~ 60 cm. Dielectric slabs line 
all four walls, their dimensions chosen to give a difference in 
insertion phase as close as possible to 90° between orthog- 
onal TEjo and TEqi modes in the waveguide. The design 
(developed with electromagnetic simulation software) gives a 
phase difference of 90° ± 2° and in practice a performance of 
90° ±5° is achieved. Th is design is based on the device devel- 
oped bv lSrikanthl d 19971) using corrugations in the waveguide 
walls to provide the differential phase shifts across the band 
18.9 to 26.5 GHz. The fractional bandwidth of the two de- 
vices is similar (essentially an entire waveguide band) and the 
performance is similar. In our case use of corrugations would 
have resulted in lower loss but the overall length would have 
been about three times larger. The dielectric that we used was 
Teflon (dielectric constant e ss 2.3). We explored designs us- 
ing materials with higher dielectric constant in the hope of 
reducing the length of the device, but were unable to make an 
improvement, either in size or performance, over the design 
using Teflon. 

Dissipative loss of the phase shifter is slightly unequal in 
the two orthogonal modes because more dielectric is used in 
one plane than in the other. The losses are about 0.05 and 
0.1 dB. 

The input and output of the power splitter are both in square 
waveguide, but the orientations differ by 45°. The device is 
fabricated from eight flat triangular plates bolted together. Its 
length is 20.6 cm. 

The orthomode transducer is a commercially available de- 
vice, originally developed for the Very Large Array by At- 
lantic Microwave (their model number OM6500). It has 
a square waveguide input and two outputs in rectangular 
WR650 waveguide. Fabricated from cast aluminum, it is a 
low-loss device (loss is probably less than 0.05 dB or 1%). 

2.3. Receiver 

The receiver is mounted at the prime focus in an insulated, 
temperature-controlled box. The air temperature inside the 
focus box is kept at 25.0 ±0.1 °C. The two-channel receiver 
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FIG. 2. — Block diagram of the receiver and signal-processing system. Par- 
allel receiver paths process L and R signals. LNA denotes low-noise amplifier 
and CAL denotes calibration noise source. 

and signal processing system are depicted schematically in 
Fig. [2] Directional couplers permit injection of noise calibra- 
tion signals. Circulators (unidirectional devices which allow 
the low-noise amplifiers to receive signals from the antenna 
ports but isolate the amplifiers from reflections from the polar- 
izer and feed) preceding the first amplifiers prevent backscat- 
tering and thus reduce cross-coupling between the two hands 
of polarization, which would lead to instrumental polariza- 
tion. The low-noise amplifiers (LNAs) are commercial un- 
cooled HEMT amplifiers with noise temperature ~35 K, but 
elements preceding them add about 30 K to their noise tem- 
perature. The receiver follows a conventional design, with a 
single-stage down-conversion taking the RF band of 1277 to 
1762 MHz to baseband. The 3-dB bandwidth of the system 
is 485 MHz. The overall system temperature, including re- 
ceiver noise, ground and atmospheric emission is estimated 
to be - 140 K. 
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TABLE 1 

Hardware specifications of the FPGA-based 
polarimeter 



ADC sampling rate 


max 2x1 GS/s 


ADC resolution 


8 bit 


Input bandwidth 


- 500 MHz 


FFT Core 


Xilinx FFT v4.0 


Number of channels (FFT) 


2 x 4096 (complex) 


Number of channels (Output) 


4 x 2048 (real) 


Precision (internal) 


2's complement, 54 bit 


Precision (output) 


2's complement, 32 bit 


Minimum integration period 


25 ms 


Telescope interface requirements 


noise source control 



3. THE FPGA-BASED POLARIMETER 

To save development time and minimize engineering effort 
the polarimeter was implemented using off-the-shelf FPGA 
and ADC cards from Sundance Multiprocessor Technology 
Ltd. 6 (UK). This approach provided an integrated hardware 
and software solution that was relatively low cost and sup- 
ported on a standard PC running Windows XP. The specifi- 
cations for the polarimeter are given in Table Q] As shown 
in Figure [5] the hardware consists of a carrier card, an ADC 
module, and two FPGA modules. The ADC module has one 
Atmel ADC and a Xilinx VP30 FPGA. The ADC has two 8- 
bit converters and samples at 1 Gsps. It outputs its data at 125 
MHz as two 8-bytes fields. These are received by the VP30 
which applies weights to the input samples in the time do- 
main (windowing) and outputs a frame of contiguous samples 
in time. These frames are distributed equally between the two 
FPGA modules on the carrier board, half going to one FPGA 
and half to the other. Each FPGA module has a Xilinx VP70, 
which first buffers the ADC fields into a frame in a FIFO (first- 
in first-out shift register). The data are clocked at a rate of 
125 MHz. Fast Fourier Transform is perfomed on frames, 
which are then summed over an integration period. The inte- 
grated power spectra are sent to the Host PC over a 33 MHz 
PCI bus. The FPGA firmware is programmed in VHDL, and 
synthesized and compiled using Xilinx ISE. The processing 
firmware is divided into separate modules and these commu- 
nicate in a standard way using FIFOs. This approach permits 
a high-level firmware management and inter-module commu- 
nication software tool, called Diamond from 3L 7 (UK), to be 
used. 

A block diagram of the polarimeter digital signal process- 
ing is shown in Figure @] The processing consists of win- 
dowing and Fast Fourier Transforming the two polarization 
streams, and then cross-multiplying and integrating these to 
produce the four polarization products (RR, LL, RL, and LR) 
as power spectra. A windowing function is applied to the data 
to minimize leakage from one channel into adjacent spectral 
channels. This is especially important for removing strong 
narrow band RFI. While a rectangular window yields the best 
frequency resolution and highest signal-to-noise ratio, spec- 
tral leakage into nearby channels is high, about -30 dB. Strong 
RFI in one channel can easily spread into many adjacent chan- 
nels and swamp the sky signal. In the current design the 
Blackman-Harris (BH) window-function is used for weight- 
ing samples in the time domain. This window function yields 

6 The system is comprised of the following components: the carrier card 
(SMT310Q), two VP70 FPGA modules (SMT398-VP70-6), a DSP module 
(SMT395-VP30-6), and an ADC module (SMT391) piggy-back on another 
FPGA module (SMT338-VP30-6). For more detailed information we refer 
to Sundance documentation or www.sundance.com 

7 www.31.com 
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FIG. 3. — A block diagram of the polarimeter. See text. 
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FIG. 4. — A block diagram of the FFT code. See text. 

a spectral leakage of less than -60 dB. This effectively sup- 
presses leakage and baseline fluctuations due to narrowband 
RFI. It is implemented as a lookup table in the VP30 on the 
ADC module. One disadvantage is that it discards 50% of the 
data and lowers the signal-to-noise ratio accordingly. In the 
future we will replace the window with p olyphase filter tech - 
niques (e.g. WOLA) as described by, e.g jKlein et al.l d2006). 

A 4096-point complex FFT is used to transform the win- 
dowed ADC data into spectra. This function is provided in 
the Xilinx supplied function library and is highly optimized. 
Because the VP70's are clocked at 125 MHz, eight FFT's 
are needed to process a 1 GHz data stream of R and L po- 
larized data. Four FFT are performed on each VP70 FPGA. 
In this scheme an FFT block receives 1/8 of the frames from 
the ADC FPGA. As well, since the R and L polarization data 
are real, one 2048 -point complex FFT can be used to compute 
the 2048-point FFT's for each polarization. The technique of 
using a complex FFT to compute two real FFTs is well known 
(e.g.lNussbaumeri l 198 lh and involves placing 2048 points of 
R polarization data in the real part of the FFT and 2048 points 
of L polarization data in the imaginary part of the FFT. After 
the FFT is complete the results for each polarization are re- 
covered yielding two 2048-point complex spectrums one for 
the R and L polarizations, respectively. 

The complex spectra are then multiplied and summed to 
yield the integrated 2048-point power spectra for the four po- 
larization products, namely RR, LL, RL, and LR. In the cur- 
rent implementation sums of 760 spectra form integrations 
with a period of 25 ms. The integrated spectra, comprising 
8192 signed four-byte integers, are sent to the host PC over 
the PCI bus. The average data rate is about 1.3 MBps which 
is easily handled on a 33MHz PCI bus. 

4. ELECTRONIC GAIN AND PHASE CORRECTION 

Amplifiers suffer from gain and phase variations which 
cause the electronic gains to vary. Phase shifts of one hand 
of polarization relative to the other can also be introduced by 
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variations in cable length due to temperature changes. These 
gain and phase drifts are corrected using a calibration signal, a 
noise signal injected into the signal path ahead of the first am- 
plifier through a directional coupler. An RF switch between 
the noise source and the coupler is used to switch the calibra- 
tion signal on and off, synchronized with the polarimeter. A 
switching time of 200 fis is is allowed for between integra- 
tions. An integration packet of 50 ms length consists of two 
integrations with cal-on and cal-off. 

The actual intensity of the calibration signal for an inte- 
gration packet is determined by calculating the difference be- 
tween the cal-on and the corresponding cal-off spectrum. The 
equivalent noise temperature of the calibration signal is ap- 
proximately 20 K. The rms noise in a single channel and a 
50 ms integration is about 8% of that and is much larger than 
the expected gain fluctuations. Therefore, the average over 
200 channels and 90s (45s of the preceding and 45s of the fol- 
lowing integrations) of the calibration signal is taken, which 
reduces the noise in the calibration signal temperature to an 
acceptable value of 3 mK. 

The calibration signal is split and equal-amplitude compo- 
nents are injected into the L and R paths. This is equivalent 
to a signal with constant polarization angle. The sky values 
are then expressed relative to the measured intensity and po- 
larization angle of the calibration signal. 

5. RFI MITIGATION 

Only a small fraction of our observing band lies within pro- 
tected allocations for radio astronomy. Hence, RFI mitiga- 
tion becomes a major, if not a limiting, factor in our obser- 
vations. RFI is present over a wide range of intensities, from 
very strong satellite emissions (GPS and GLONASS among 
others), to extremely weak RFI which is seen only in the final 
data processing steps after sufficient averaging. The same can 
be said about the time-domain: in some channels RFI can be 
present at all times, while it is sporadic in others. In order 
to minimize the level of locally generated RFI the observa- 
tory buildings are screened, which attenuates radio signals in 
L-band by about 25 dB. Known sources of RFI at the Obser- 
vatory, such as computers, printers, laboratory equipment, and 
digital signal-processing equipment are enclosed in shielded 
boxes or are operated inside Faraday-cage screened rooms to 
further suppress their emission (by at least 50 dB). 

We found that RFI is more severe at low elevations when 
the telescope points south than at the Zenith. This may be 
due to two effects: geostationary satellites lie at low elevation 
as seen from Canada, and the spillover sidelobes look at the 
observatory building and pick up terrestrial RFI, most likely 
produced by the surrounding community. In contrast, with 
the telescope at the Zenith the spillover lobes see the ground. 
Rather than to try to correct for RFI in any way, our approach 
is to flag any RFI found in the data and use the remaining 
(interference-free) part of the spectrum for data analysis. We 
applied a three-stage RFI search to the data at different places 
in the data processing chain, reflecting the increase in sensi- 
tivity due to averaging. The first two stages are part of the 
real-time processing chain, and the last stage is applied to the 
final maps. If RFI is detected in one or more correlation prod- 
ucts, all four correlation products are flagged for that channel. 

The first algorithm detects RFI in the time domain. It de- 
tects strong, time-variable signals and flags affected channels 
until 60 seconds after the last detection. On average about 
10 to 20 channels are flagged in this process. A high thresh- 
old at this stage is necessary to avoid erroneous flagging of 



o.i 



o.oi 




10 100 1000 

Integration Time (s) 



FIG. 5. — Allan-variance plot of total power data in a single channel at 
1460 MHz, measured with the telescope pointing towards the North celestial 
pole (continuous line). For reference the Allan variance of white noise is 
plotted (dashed line). The receiver is stable over periods of at least 500s. 
Longer term drifts are probably caused by changes in the level of Galactic 
emission received through the sidelobes during the measurement. 

good data (for example, while scanning across a bright point 
source which raises the baseline rapidly). The second algo- 
rithm works in the frequency domain. A simple median fil- 
ter is used to detect signals above or below a certain thresh- 
old. Naturally this algorithm is more successful than the time- 
domain search and flags about 150 to 250 channels (about 
10% of the band) at any time. A third RFI search is applied 
to the final maps using another median filter algorithm. It is 
usual for an additional 10% of the band to be flagged during 
this stage, which leaves about 80% of the band usable for ra- 
dio astronomy. 

Fortunately, the Observatory site is in mountainous coun- 
try and terrain shielding very effectively attenuates interfering 
signals from terrestrial sources. Interference from satellites 
and airborne transmitters is not overwhelming. With the im- 
plementation of the three relatively simple RFI excision pro- 
cedures it has been possible to make sensitive observations 
over a band in which radio astronomy has little protection 
from other spectrum users. 

6. STABILITY & POLARIZATION PERFORMANCE 

The short-term stability of the whole system was evaluated 
by an Allan- Variance measurement with the telescope point- 
ing towards the North celestial pole. The Allan variance plot 
(Fig. |5]l shows that received noise behaves like white noise 
and no systematic drifts occur on time scales of less than 500 
seconds. Variations on timescales longer than 500 seconds 
probably arise from changing signal levels as the sky drifts 
through the telescope sidelobes. 

The long-term stability of the system was examined based 
on the daily calibration sources. A variation in gain of the 
whole system of 1 % over periods of one month was found, 
possibly caused by small temperature variations inside the fo- 
cus box. These small instabilities are removed by calibration 
against astronomical sources. 

Instrumental polarization has two components, cross-talk 
between R and L channels (occuring in the feed and waveg- 
uide components and their interaction with reflections at the 
low-noise amplifier inputs). Cross-talk varies from 25% at the 
bottom end of the band to < 3% towards the high frequency 
end. The cross-polar response of the antenna has not been 
taken into account in the results presented here; this effect 
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is evident only where total intensity is very high - along the 
Galactic plane. 

7. PROOF-OF-CONCEPT STUDY 

The first 8% of data (300 of a total of 3600 scans) of the 
GMIMS component survey with the DRAO 26-m Telescope 
are used here for a proof-of-concept study. The goal of this 
study is to develop processing, calibration, and analysis al- 
gorithms, which will eventually be applied to the complete 
DRAO survey as well as to the other component surveys car- 
ried out by the GMIMS collaboration. The result of this study 
is the first RM-Synthesis of Galactic, diffuse polarized emis- 
sion, obtained with a single-antenna telescope. 

7.1. Observations 

A set of scans along the Meridian, 12 arcmin spaced, was 
defined which would provide sampling slighly surpassing the 
Nyquist criterion of the entire sky from —30° to 87° decli- 
nation at the highest frequency. Long scans extend from the 
horizon limit at —30° declination to 87°, close to the North 
celestial pole, and short scans from —30° to 60°. To avoid 
oversampling around the pole every second scan is a short 
scan. A quasi-random sequence of scans was devised, aiming 
to avoid systematic effects while making optimal use of the 
observing time. With a scanning speed of 52.5 arcmin/min 
short scans take 1.7 hours and long scans 2.2 hours. Scans 
are made exclusively at night to avoid contributions from so- 
lar emission through the sidelobes, but calibration sources are 
observed during the hour before sunset and after sunrise. For 
the final survey each scan will be observed twice, as an up and 
a down scan, resulting in double coverage of every pixel. The 
region presented in this paper was observed in April 2008 and 
is non-uniformly covered. 

Resolution was reduced to 1° by smoothing in the image 
plane to approximate full Nyquist sampling. The rms noise 
in one channel (237 kHz width) was found to be ~ 80 mK in 
total power and ~ 40 mK in polarization. This is about 1 .6 
times higher than the theoretical value of 25 mK, calculated 
based on our system temperature, bandwidth, and integration 
time. The reasons for the higher rms noise are a) the deac- 
tivation of one of the four FFT lanes on each FPGA due to 
temperature problems 8 , which resulted in a loss of 25% of 
digitized samples and thus increased the rms noise by a fac- 
tor of 1.15; and b) weighting of the digitized samples using 
a Blackman-Harris windows function, which, in effect, made 
use of only 50% of the input signal and thus increased rms 
noise by a factor of 1.41. 

7.2. Post-Processing & Calibration 

We now discuss the most important post-processing and 
calibration steps. These are: 1) correction of instrumental ef- 
fects (including corrections for bandpass, relative phase shift 
across the band, and cross-talk introduced in the receiving 
system by a left polarized signal entering the right channel 
and vice-versa); 2) subtraction of ground radiation; and 3) re- 
moval of scanning effects. A correction of the cross-polar re- 
sponse of the telescope beam, which requires a deconvolution 
of the final maps with the beam pattern, was not attempted. 
While bandpass and phase corrections were averaged over a 
few channels, all other corrections were performed on each 
channel individually. 

8 The problem has since been fixed and all four lanes are now in operation. 
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FIG. 6. — Relative gain of RR (continuous) and LL (dashed) on the left 
axis, as obtained through a single observation of Tau-A. The right axis shows 
the phase difference (RR — LL) across the band as measured through a series 
of 3C286 and 3C270 observations. Plotted are the actual measurement (gray) 
and a smoothed version (thin black) that is used for the calibration. The phase 
is set to zero at the highest frequency. 



TABLE 2 

Primary and secondary gain calibrators 



Name 


Literature (1.4 GHz) 


Adopted (1.4 GHz) 


Remarks 




Flux a 


Flux a 




Cyg-A 


1579 Jy -1.01 


1589 Jy -1.07 


1,3 


Tau-A 


908 Jy -0.28 


848 Jy -0.27 


1.3 


Vir-A 


208 Jy -0.83 


207 Jy -0.90 


1.3 


Cas-A 


2442 Jy -0.78 


1861 Jy -0.77 


2,4 



"Primary Calibrator. 
^Secondary Calibrator. 

c Flux value taken from the VLSS Bright Source Spectral Calcula- 
tor (Helmboldt et al. 2008), which is based on Baars et al. (1977). 

rf The literature flux value is the value for the epoch 1980 from 
Baars et al. (1977). 

In addition to the electronic gain calibration, which is part 
of the real-time processing chain (Sect. |4|, daily flux calibra- 
tions are derived from observations of two of our four flux 
calibrators (see Fig.|6]). We used Cyg-A, Tau-A, and Vir-A as 
our primary calibrators. Cas-A was used as secondary cali- 
brator whenever none of the primary calibrators was visible. 
The flux of the secondary calibrator was derived from the pri- 
mary calibrators. A calibration observation consisted of 1 1 
declination scans spaced by 12 arcmin, yielding a map of the 
source. A two-dimensional Gaussian was fitted to the cali- 
bration observations and the amplitude of the fitted Gaussian 
was used as the calibration value. Table|2]summarizes the flux 
values and spectral indices (a) used for the primary calibra- 
tor and derived for the secondary calibrators. We found that 
flux values slightly different from the published values gave 
the best overall consistency across the four calibrators. The 
final data are tentatively calibrated in main-beam brightness 
temperature by assuming an aperture efficiency of 53% over 
the entire band, a value originally determined for 1 .4 GHz by 
iHiggs & Ta pping (2000). In the future, when survey data over 
a much larger period need to be calibrated, it is planned to take 
intrinsic variability of Cas-A into account. 

The relative phase between the two polarization channels 
(R and L) was unknown and was determined through observa- 
tions of 3C286 and 3C270. Both sources are at high Galactic 
latitudes and have well known RMs of w 0. At 1 .4 GHz, the 
percentage polarization and polarization angles of 3C286 and 



7 



3C270 are 9.3% and 3 2°, and 8.0% and 123°, respectively 
dTabara & Inoud[T980h . This measurement revealed a phase 
gradient of about 30° across the band (Fig. [5). A phase cali- 
bration is not done on a regular basis since this gradient will 
be stable until changes to the receiver are made that affect the 
path length or delay of the signal paths. 

Cross-coupling of one hand of polarization into the other 
results in instrumental polarization. In the final polarization 
maps this effect leads to a conversion of Stokes I into U and Q, 
leading to the additional signals Ui ns t = /u I and Qj nst = /q I. 
This type of instrumental polarization is the result of imper- 
fections in the feed and mismatches in the receiver. The coef- 
ficients f\j and /q were determined daily from the flux calibra- 
tion. The integrated polarization from the primary calibrators 
is very close to zero and they were considered to be unpolar- 
ized. Any polarization detected is then of instrumental origin. 
As it is evident in the final maps, the removal of cross-talk is 
not sufficient for a complete correction for instrumental polar- 
ization. The images still show spurious polarization features 
around bright sources produced by the cross-polar response 
of the telescope. The correction of cross-polar response will 
require a cleaning (de-convolution) of the final polarization 
data, which is planned for the complete survey and requires 
fully Nyquist-sampled maps in all Stokes parameters. 

Ground radiation is received through the side- and spillover 
lobes of a telescope. It can contribute a signal an order of 
magnitude stronger than the sky, especially at low elevations. 
As a temporary approach for the data presented in this paper, 
ground radiation profiles were determined by plotting the ob- 
served Stokes I, U, and Q intensities versus elevation. Only 
scans made between 13h and 17h R.A. were used, avoiding 
the bright emission from the Galactic plane. Ground contri- 
bution was then retrieved from these scatter plots by taking 
the lower envelope in Stokes I, and a median value for Stokes 
U and Q. 

Scanning effects (stripes) in the final maps were removed 
by a technique known as "basket weaving". This technique 
removes baseline effects by comparing the overlaps of scans 
made in different directions and times with the assumtion 
that th e baselines vary a nd the sky brightness remains con- 
stant (ISieber et al.lll979l) . The stripes are most likely caused 
by variations of the electromagnetic properties of the ground 
surrounding the telescope affecting the ground pick-up, time- 
variability of the noise figure of the LNAs, or small errors 
in the flux calibration. Basket weaving was applied on each 
channel individually. To cope with the high data volume the 
basket weaving algorithm ran on a cluster consisting of 40 
CPUs. A problem with the basket weaving technique is the 
floating baseline level which can introduce large-scale distor- 
tions in the maps. The Stokes I map is more affected by these 
distortions then the Stokes U and Q maps. The maps pre- 
sented here (Fig. |7]i show the R.A. range mostly free of dis- 
tortions (a gradient is still present in the total intensity map at 
R.A. > 20h). 

7.3. Application of Rotation Measure Synthesis 

RM-Synthesis consists of two separate steps: 1) de-rotation 
of Stokes U and Q, and 2) cleaning of the RM-Spectra. In the 
first step, the final data cube is de-rotated for a range of RM 
values. Here we used RMs from -1250 to 1250 rad m~ 2 in 
steps of 10 rad vcT 2 . For each RM value the Stokes U and 
Q values in each pixel and channel are de-rotated by the an- 
gle RMA 2 . The de-rotated cube is then collapsed (averaged) 
into a 2-D map of polarized intensities, which forms a slice of 
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FIG. 7. — Maps show the first 8% of survey data, tentatively calibrated as 
explained in the text. Total intensity is shown at 1300 MHz, and polarized 
intensity is shown at 1300 and 1700 MHz. Contour lines in total intensity 
correspond to 500, 1000, 1500, 2000, 4000, and 8000 mK. In polarized in- 
tensity contour lines indicate polarized intensity in steps of 100 mK at 1300 
MHz, and 50 mK at 1700 MHz, resp. Features in the maps are explained in 
the text. 

the RM-Synthesis cube. RM-Spectra are constructed for each 
pixel from the series of 2-D maps at different RMs. 

The second part of RM-Synthesis is cleaning. Cleaning is 
necessary because the observed spectra are "contaminated" 
with spurious line features caused by sidelobes in the rotation 
measure spread function (RMSF). The RMSF is the response 
function (the instrumental response) with which the spectrum 
is convolved and is the result of incomplete frequency cover- 
age. Figure[8]shows an RMSF for our observations, which is a 
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FIG. 8. — A cleaned (solid line) and uncleaned (dotted line) RM-Spectrum 
for a position close to the centre of the surveyed region. The uncleaned spec- 
trum represents a typical RMSF for our data. 



TABLE 3 
RM-Synthesis Parameters 



Parameter 


Value 


Resolution in Faraday Depth 


132 rad m" 2 


Largest Scale in Faraday Depth 


108 rad m -2 


Maximum Observable Faraday Depth 


10 5 rad m- 2 



typical RMSF for data in our frequency range, coverage, and 
flagging. Using an analogy to synthesis imaging, the RMSF 
can be understood as the "dirty beam". The observed spectra 
("dirty images") need to be de-convolved (with an appropri- 
ate clean algorithm) to obtain the "cleaned images". The clean 
model consists of a series of delta peaks. To obtain cleaned 
RM-spectra the RM-model is convolved with a Gaussian (the 
"restoring beam") with a width corresponding to the expected 
resolution in RM-Synthesis, which, in our case, was slightly 
less than 132 rad m~ 2 . Figure [8] shows a cleaned and an un- 
cleaned RM-Spectrum. Spurious components at the ^5% 
level are still present and are visible at ±850 rad mT 2 in the 
spectrum shown. 

Our frequency range and spectral resolution determine: 1) 
the resolution of the RM-Spectra (the resolution in Faraday 
depth), 2) the largest scale in Faraday depth we are sensitive 
to, and 3) the maximum observable Faraday depth. These val- 
ues are summarized in Tab. [3] In our case the largest scale in 
Faraday depth we can measure is smaller than the resolution. 
This is an important constraint to keep in mind: two peaks 
in our RM-spectrum do not necessarily indicate two different 
Faraday rotation layers, but may also arise from a continuous 
distribution of RMs in Faraday depth, a few RMTFs wide, of 
which only the two edges are resolved. 

7.4. General Discussion 

Figure [9] shows two slices of the RM-Synthesis cube. The 
rms sensitivity of these images is ~ 3 mK, the rms noise ex- 
pected from the entire 470 MHz bandwidth. This high sen- 
sitivity is one of the great benefits of RM-Synthesis. Fur- 
thermore, the abundant frequency coverage permits quite ag- 
gressive flagging of RFI, leading to good signal-to-noise ratio 
even outside the protected radio astronomy bands. Finally, we 
note that RM-Synthesis eliminates the problem of bandwidth 
depolarization: it provides an adaptive filter that can be tuned 



rad m 2 




R.A 

+60 rad m 2 




20h lgh 16h 14h 12h 

R.A. 



FIG. 9. — The upper map shows the RM-Synthesis slice at rad m 2 with 
contour lines every 100 mK of polarized intensity. The lower map shows the 
60 rad m — 2 slice with contour lines in steps of 25 mK. 

to the characteristics of the emission. 

The strength of RM-synthesis, however, is in revealing 3- 
D structures in Faraday space. This requires a A 2 -coverage 
matched to the expected Faraday rotation. Our resolution in 
Faraday depth is probably sufficient to map the general trend 
of RMs, but will only allow us to resolve out multiple RM 
components for highly Faraday rotated parts of the ISM. Our 
frequency range (A 2 -coverage), however, will allow us to de- 
tect highly Faraday rotated structures even if the fractional 
polarization of these regions is extremely low. Without RM- 
Synthesis the data would be "confusion limited" in the sense 
that faint but highly Faraday rotated emission could not be 
separated from the usually stronger, unrotated "foreground". 
RM-Synthesis helps by separating differently Faraday rotated 
structures. 

The data presented here are preliminary products in the 
sense that they are not yet absolutely calibrated: ground ra- 
diation profiles were determined in an area of the sky contain- 
ing large-scale Galactic emission in Stokes I, U, and Q, which 
makes the separation of ground and sky emission difficult. An 
independent, absolute temperature scale has not been estab- 
lished for the telescope. Instrumental polarization is visible in 
the maps caused by the cross-polar response of the telescope, 
which is most obvious along the Galactic plane. Smaller prob- 
lems are remaining scanning effects. Despite the preliminary 
state of these data, the final Stokes U and Q maps at 1 .4 GHz 
compare very well with t he DRAO Low R esolution Polariza- 
tion Survey at 1.4 GHz (IWolleben et al.ll2006l) . which con- 
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vinces us that our proof-of-concept study will provide valu- 
able input for future RM-Synthesis observations with single- 
antenna telescopes. 

7.5. Polarized Emission at rad mT 1 

A few large-scale structures in the rad m~ 2 RM-Synthesis 
slice (Fig.|9]top) are immediately obvious. Polarized emission 
from the North-Polar Spur is clearly visible as an arc-like fil- 
ament in the centre of the map, going in R.A. from 14h to 
17h30m, and in DEC from -5° to 20°. South of the Galac- 
tic plane, at R.A. from about 19h to 21h30m and DEC from 
—25° to 0°, a patch of polarized emission is visi ble, which 
has b een interpreted as part of the Loop I structure ( Wolleben 
120071) . At R.A. = 18h30m and DEC = 60° another polarized 
filament is visible, most likely associated with Loop III. Polar- 
ized emission within a narrow strip along the Galactic plane 
can be assumed to be instrumental polarization. 

7.6. Polarized Emission at High Positive RMs 

The RM-Synthesis slice at 60 rad m -2 (Fig. [9] bottom) is 
chosen as an example for polarized emission at high positive 
RMs. Most features in this slice actually peak at slightly dif- 
ferent RMs. The most obvious structure is centred on R.A. = 
17h30m and DEC= 20°, about 20° in size. Other structures 
along the Galactic plane and structures close to the edges of 
the map may be spurious. 

Previous polarization surveys could detect polarized emis- 
sion from the North-Polar Spur only at Galactic latitudes 
above ~ 25°, while in total intensity the spur can clearly be 
traced down to lower latitudes. The two filaments found at 
high positive RMs run from the southern tip of the polarized 
emission of the North-Polar Spur to the Galactic plane. The 
two filaments may therefore be associated with the North- 
Polar Spur, perhaps representing highly Faraday-rotated low- 
latitude emission of the spur. We postpone further analysis 
and interpretation until sampling and angular resolution are 
improved. 

In summary, the two filaments at 60 rad m~ 2 are good ex- 
amples of structures that are only visible in polarization and 
which can be revealed only by RM-Synthesis. Such high RM 
require a strong magnetic field parallel to the line-of-sight. 
The filaments are most likely structures in the magnetic field, 
perhaps produced by supernovae or stellar winds in the inter- 
stellar medium, compressing the ambient magnetic field. 

7.7. Polarized Emission at Negative RMs 

Our RM-Synthesis data reveal faint, Faraday-rotated struc- 
tures at RM < —25 rad m~ 2 . Some of these structures are 
filamentary in shape and seem to be associated with the North- 
Polar Spur, although not with its southern tip but other parts 
of its shell at higher Galactic latitudes. Thus, it is possible that 
our RM-Synthesis data reveal, in fact, the 3-D structure of the 
magnetic field of the North-Polar Spur, but more complete 
coverage of this region is required to confirm this interpreta- 
tion. 

7.8. Conclusions from the Proof-of-Concept Study 

Wide-band polarimetric data present a challenge for cali- 
bration in that the receiver characteristics and performance 
can vary significantly over the band. Furthermore, the orien- 
tation and magnitude of sidelobes and spillover lobes depends 
on frequency, resulting in highly frequency-dependent ground 



radiation profiles. In our data we find channels that are virtu- 
ally free of instrumental polarization while other channels are 
almost free of ground radiation. 

These effects can be very useful for calibration. For ex- 
ample, channels free of ground radiation can be used as "an- 
chor points" for an absolute calibration, making the removal 
of ground radiation more robust in the presence of large-scale 
emission structures in the data. Similarly, channels free (or 
nearly so) of cross-talk can be very useful for the determina- 
tion of the cross-polar component of the instrumental polar- 
ization. 

An important observation is that we find unexpectedly high 
RM values for the diffuse emission. Diffuse emission is sub- 
ject to various depolarization effects along the line-of-sight, 
which has led to the argument that most of the diffuse polar- 
ized emission at < 2 GHz must be of local origin within a few 
kpc. We find RMs as high as 100 rad m~ 2 for some of the 
diffuse emission, which may either indicate that the emission 
comes from a larger distance, or that the local magnetic field 
is stronger than previously believed. 

Previous rotation measure surveys of the diffuse Galactic 
emission revealed only moderate RMs of less than a few 
tens of rad m~ 2 dBrouw & Spoelstral [l976). We conclude 
that Faraday rotation should be measured using wide-band 
spectro-polarimetric observations. The relation between po- 
larization angle and A 2 is not always linear for the diffuse 
Galactic emission, rendering the traditional techniques of de- 
termining RM by least-square fitting inadequate for diffuse 
emission. 

8. SUMMARY & OUTLOOK 

Recent developments in digital signal processing allowed 
us to build a digital polarimeter in a short time. In this pa- 
per we described the technical process from building a digi- 
tal polarimeter, designing the feed, and obtaining, calibrating, 
and analyzing wide-band polarization data. For the first time, 
RM-Synthesis was successfully applied to wide-band spectro- 
polarimetric observations of the diffuse Galactic emission ob- 
tained with a single-antenna telescope over a large region on 
the sky. 

Observations for the DRAO Rotation Measure Survey are 
ongoing and are 42% completed as of September 2009. The 
final survey will be better than Nyquist sampled with an rms 
sensitivity, in a 12 arcmin wide pixel, of 105 mK per channel 
or 3 mK for the entire band (taking into account data loss due 
to RFI). Although not yet in a publishable state, a first look 
at these data confirms the presence of the polarized structures 
discussed in this paper. 

This survey forms part of GMIMS. Three of the six compo- 
nent surveys of the GMIMS project are now underway, with 
another two surveys in the planning stages. An important ob- 
jective of GMIMS is that all measurements be absolutely cali- 
brated. In practical terms this means that the calibration noise 
signal will be carefully measured relative to terminations at 
cryogenic temperatures (liquid nitrogen), and that the aperture 
efficiency will be determined using calibrators of well-known 
flux density. 

The authors would like to thank Sundance and 3L for out- 
standing support during the development phase of the digital 
polarimeter. We are grateful to the IT staff at DRAO, Tony 
Hoffman and Peter Cimbaro, for providing and maintaining 
the telescope infrastructure. We also thank Ron Casorso, Paul 
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scope. The Dominion Radio Astrophysical Observatory is a Canada. 



REFERENCES 



Baars, J.W.M., Genzel, R., Pauliny-Toth, I.I.K., & Witzel, A., 1977, A&A, 

61,99 

Bernardi, G., et al. 2009, arXiv:0904.0404 
Brentjens, M.A., & de Bruyn, A.G. 2005, A&A, 441, 1217 
Brouw, W. N., & Spoelstra, T. A. T. 1976, A&AS, 26, 129 
de Bruyn, A.G, & Brentjens, M.A., 2005, A&A, 441, 931 
Burn, B. J. 1966, MNRAS, 133, 67 

Carretti, E., Haverkom, M., McConnell, D., Bernardi, G, Cortiglioni, S., 
McClure-Griffiths, N. M., & Poppi, S. 2008, arXiv:0806.0572 

Gibson, S. J. 2005, Bulletin of the American Astronomical Society, 37, 1489 

Helmboldt, J. R, Kassim, N. E., Cohen, A. S., Lane, W. M., & Lazio, T. J. 
2008, ApJS, 174,313 

Higgs, L. A., & Tapping, K. F. 2000, AJ, 120, 2471 

Klein, B., Philipp, S.D., Glisten, R., Kramer, I., & Samtleben, D. 2006, 
Proc. SPIE, 6275 

Pacholczyk, A. G. 1970, Radio Astrophysics, Series of Books in Astronomy 

and Astrophysics, San Francisco: Freeman, 1970 
Nussbaumer, H. J. 1981, Fast Fourier Transform and Convolution 

Algorithms, Springer- Verlag Berlin Heidelberg New York (Springer Series 

in Information Sciences) 



Reich, W., Fiirst, E., Reich, P., Uyaniker, B., Wielebinski, R., & Wolleben, R. 

2004, The Magnetized Interstellar Medium, 45 
Schnitzeler, D.H.F.M., Katgert, P., Haverkom, M., & de Bruyn, A.G, 2009, 

A&A, 461, 963 

Schnitzeler, D.H.F.M., Katgert, P., & de Bruyn, A.G, 2009, A&A, 471, L21 
Schnitzeler, D.H.F.M., Katgert, P., & de Bruyn, A.G., 2009, A&A, 494, 611 
Sieber, W., Haslam, C. G. T., & Salter, C. J. 1979, A&A, 74, 361 
Smegal, R. J., DRAO Technical Memo, 2005 

Srikanth, S. 1997, IEEE Microwave and Guided Wave Letters, 7, 150 
Tabara, H., & Inoue, M. 1980, A&AS, 39, 379 
Testori, J. C, Reich, P., & Reich, W. 2008, A&A, 484, 733 
Wohlleben, R., Mattes, H., & Lochner, O. 1972, Electronics Letters, 8, 474] 
Wolleben, M., Landecker, T.L., Reich, W., & Wielebinski, R. 2006, A&A, 
448,411 

Wolleben, M. 2007, ApJ, 664, 349 

Wolleben, M., et al. 2008, arXiv:0812.2450 



